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An analytical variational method is applied to the molecular Holstein Hamiltonian in which the 
dispersive features of the dimension dependent phonon spectrum are taken into account by a force 
constant approach. The crossover between a large and a small size polaron is monitored, in one, two 
and three dimensions and for different values of the adiabatic parameter, through the behavior of 
the effective mass as a function of the electron-phonon coupling. By increasing the strength of the 
inter-molecular forces the crossover becomes smoother and occurs at higher e-ph couplings. These 
effects are more evident in three dimensions. We show that our Modified Lang-Firsov method starts 
to capture the occurence of a polaron self-trapping transition when the electron energies become of 
order of the phonon energies. The self-trapping event persists in the fully adiabatic regime. At the 
crossover we estimate polaron effective masses of order ~ 5 — 40 times the bare band mass according 
to dimensionality and value of the adiabatic parameter. Modified Lang-Firsov polaron masses are 
substantially reduced in two and three dimensions. There is no self-trapping in the antiadiabatic 
regime. 

PACS: 71.38.+i, 63.10.+a, 31.7Q.Ks 



I. INTRODUCTION 



There has been a growing interest towards polarons 
over the last years also in view of the technological po- 
tential of polymers and organic molecules [1,2] in which 
polaronic properties have been envisaged. Theoretical 
investigations on polarons generally start from the Hol- 
stein Hamiltonian [3] originally proposed for a diatomic 
molecular chain along which hopping of electrons, lin- 
early coupled to the vibrational quanta, takes place ac- 
cording to a tight-binding description. If the local e-ph 
coupling is sufficiently strong the induced lattice defor- 
mation may dress the electron and transform it into a 
polaronic charge carrier [4,5]. The conditions for polaron 
formation and its mobility properties may however de- 
pend also on the adiabaticity ratio, on dimensionality, on 
peculiarities and anharmonicities of the lattice structure 
[6-15]. As the Holstein Hamiltonian can be identically 
applied to states made of excitons [16,17] and phonons 
it also provides a useful tool in optical spectra analysis 
and transport properties of organic materials [18,19] with 
large scale applications. 

While the physical properties of polaronic systems 
change [20,21] according to the size of the polaronic 
quasiparticle, a number of theoretical tools [22-26] has 
been applied to clarify nature and width of the crossover 
between a large (with respect to the lattice constant) 
polaron at weak e-ph coupling and a small polaron at 
strong coupling for a given value of the adiabatic pa- 
rameter. Provided that a phase transition is ruled out in 
the Holstein Hamiltonian with dispersive optical phonons 
[27] being the ground state energy analytic in the e-ph 
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coupling, such a crossover may still appear as a smooth 
transition in the antiadiabatic regime or rather as a sud- 
den (but continuous) event in the adiabatic regime [28]. 
While the narrowing of the polaron band signals the on- 
set of the crossover it is certainly from the behavior of 
the effective mass that the self-trapping event can be ac- 
curately located in the intermediate region of e-ph cou- 
plings [29]. Although precisely in the latter region per- 
turbative studies traditionally reveal their shortcomings 
an analytical method, the Modified Lang-Firsov (MLF) 
transformation, [30] has been developed to overcome the 
limitations of the standard Lang-Firsov (LF) [31] ap- 
proach on which strong coupling perturbation theory is 
based. As an enhancement of the polaron mass should 
be accompanied by a reduction in the polaron size, the 
electron-phonon correlation function offers an indepen- 
dent tool to analyse the crossover through the measure 
of the spread of the lattice deformation. Since the no- 
tion of self-trapping transition has often assumed differ- 
ent meanings in the literature we emphasize that our 
view of a self-trapped polaron is not that of a localized 
and immobile object but, rather, of a small quasiparticle 
whose different ground state properties have undergone 
a transition, driven by the e-ph coupling, at distinct al- 
though closely related points in the polaron parameter 
space [32]. Among these properties we study in this pa- 
per, using the MLF transformation [33,34], the polaron 
mass and the correlation function as obtained from a Hol- 
stein Hamiltonian in which the dispersion of the optical 
phonon branches is fully accounted for in any dimen- 
sionality. Besides depicting a model more appropriate 
to physical systems, dispersive phonons represent a rel- 
evant feature of the Holstein model itself as previously 
shown by one of us [35]. The role of the intermolecular 
forces in the crossover of the MLF polarons at different 
dimensionalities is a main focus of our investigation. The 
generalities of the model are given in Section II while the 
results are presented in Section III both for the polaron 
mass and for the static correlation function. Section IV 
contains some final remarks. 



II. MODIFIED LANG-FIRSOV PHONON BASIS 
FOR THE HOLSTEIN MODEL 

We consider the dimension dependent Holstein Hamil- 
tonian consisting of one electron hopping term, an inter- 
action which couples the electronic density and the ionic 
displacements at a given site and dispersive harmonic op- 
tical phonons as: 

^ = E + 9 E "^(4 + + E'^q^q^q (0.1) 

<ij> i q 

the first sum is over z nearest neighbors, cj and Cj are 
the real space electron creation and annihilation opera- 
tors. Hi (= ctcj) is the number operator, &t and 6j are 
the phonons creation and annihilation operators. 6^ is 
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the Fourier transform of bl and is the frequency of 
the phonon with vector momentum q. 

The standard practice in dcahng with the Hamiltonian 
(1) is to apply the LF transformation where a phonon ba- 
sis of fixed displacements (at the electron residing site) 
is used. Such a choice of phonon basis diagonalizes the 
Hamiltonian in absence of hopping. The hopping term 
is then treated as a perturbation. However the LF ap- 
proach under simple approximations, e.g. within zero- 
phonon averaging or zeroth order of perturbation, can- 
not describe the retardation between the electron and 
the lattice deformations produced by the electron. This 
retardation induces a spread in the size of the polaron 
and becomes very important for weak and intermediate 
e-ph coupling. The MLF phonon basis, where the dis- 
placements of the oscillators at different sites around an 
electron are treated variationally, can describe the retar- 
dation and a large to small polaron crossover even within 
simple approximations [33,34]. Recently the convergence 
of the perturbation series within the LF and the MLF 
methods has been studied in a two-site Holstein model 
for the ground state [30] as well as for the first excited 
state [36]. It was found that: (i) within the MLF method 
the perturbation corrections are much smaller than those 
corresponding to the LF method in the range from weak 
to intermediate e-ph coupling, (ii) the convergence of the 
perturbation series within the MLF is also much better 
in that range, (iii) in the strong coupling limit the MLF 
phonon basis reduces to the LF basis and the LF pertu- 
bation method works very well in this limit. The above 
studies have clearly pointed out that the MLF pertur- 
bation method works much better than the LF method 
when the entire range of the e-ph coupling is considered. 

The MLF perturbation method has also been ap- 
plied to a many-site Holstein model with dispersionless 
phonons in ID and the supremacy of the MLF method 
over the LF method in predicting the ground state energy 
and dispersion of the polaron has been observed [37] . For 
the present case of dispersive phonon we apply the MLF 
transformation to the dimension dependent Hamiltonian 
(1): 

H = e^He~^ (0.2) 

where 

q 

= E -^^""'^^ = ^T. (0.3) 

and AqS are the variational parameters which represent 
the shifts of the equilibrium positions of the oscillators 
(quantized ion vibrations) with momentum q. For con- 
ventional Lang-Firsov transformation Aq = g/cOq- The 
MLF transformed Hamiltonian for a single electron case 
is obtained as 
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H = -CpY^m-tpY^ c\cj 

i ij 

^ q 

^ q 

+ ^q^q'^'i + - Kuj^)n^{bl^ + 6q) (0.4) 



q 

where 



= ^ 51(25' - Aqa;q)Aq (0.5) 



q 



is the polaron self-energy and 

t,=iexp[-l^A^(l-^)] (0.6) 



TV 

q 



is the polaronic hopping. The coordination number z 
is twice the system dimensionality. 

where i and j are nearest neighbor sites. As unperturbed 
Hamiltonian we choose Hq as 

Ho = -ep ^ ni + ^ ^cp\b^ (0.7) 



The remaining part of the Hamiltonian (H—Hq) in the 
MLF basis is considered as the perturbation part. The 
energy eigenstates of Hq are given by 

l'/'i-{"q}> = cl\0)e\nq,,nq^,nq^,...)ph (0.8) 

where, i is the electron site and riq^ , , are the 
phonon occupation numbers in the phonon momentum 
states qi,q2,q3, respectively. The lowest energy cigen- 
state of the unperturbed Hamiltonian has no phonon ex- 
citations, i.e. riq = for all q. The ground state has an 
energy — —ep and is A^-fold degenerate, where is 
the number of sites in the system. The perturbation lifts 
the degeneracy and to first order in t the ground state 
energy of the 3D- polaron with momentum k is given by 

£;o(k) = -Cp - tpju (0.9) 

and the corresponding eigenstate is |k, nq = 0) = 

The second order correction to the ground-state energy 
of the polaron with momentum k is given by 



= E E ^ ^ ^ , |({nq},k'|g - Ho\k, {0})^ 

(0.10) 
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It is evident that the second order correction has con- 
tributions from intermediate states having all possible 
phonon numbers, i.e. each in Eq. (10) takes values 
from zero to infinity with the condition that titot = 

En > 1. 

By minimizing the zone center ground state energy we 
get the variational parameters Aq: 

\ = i- (0.11) 

and, by such a choice of Aq, the one phonon matrix 
element between the ground state |k = 0, {riq = 0}) and 
the first excited state 



(lq,k'|#-ffo|k = 0,{0}) =5k^ 

'~y 

[-ZtpXq{l - y ) + (5 - AqWq; 



1 

X 

N 



(0.12) 



vanishes. Then, the one phonon excitation process 
yields no contribution to the second order correction for 
the MLF ground state energy. The Aq's appropriate to 
the ID, 2D and 3D systems are easily obtained by (11). 

The static correlation function involving the electron 
charge at i-th site and the lattice deformation at the i+n 
site are given by 

Xn = {i>G\cUi{bl+n + bi+n)\i^G)/2g{ni) (0.13) 



where g = N~^^^{g/ojq,) and \iJjg) denotes the 
ground state for the polaron with momentum k=0. The 
denominator in Eq. (13) is used to normalize the cor- 
relation fimction with respect to its on-site value in the 
strong coupling limit. Ui is the electron number operator 
and {ni)= 1/N for the 1-electron system. While Cj and 
bi are the bare electron and phonon annihilation opera- 
tors in the undisplaced oscillator basis respectively, the 
corresponding operators in the MLF basis are the anni- 
hilation operators for the polaron and that of the phonon 
in the variationally displaced oscillator basis. The cor- 
relation function is calculated in the MLF basis within 
zero phonon averaging. 



III. POLARON MASS AND CORRELATION 
FUNCTIONS 

Previous investigations have pointed out that the Hol- 
stcin model with a dispersionlcss spectrum {loi = 0) or 
with weak intermolecular forces (wi <C u)o) would pre- 
dict larger polaron bandwidths in lower dimensionality 
against physical expectations [35]. Moreover, as pointed 
out by Holstein in his original papers [3], dispersionless 
phonons would lead to an unphysical divergent site jump 
probability for the polaronic quasiparticle [38]. Hence, 
intermolecular forces are a key ingredient of the Holstein 
model. 
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Numerical analysis [35] has shown that the bandwidths 
Ai?rf grow faster versus the intermolccular energy lui in 
higher dimensionality d thus providing a criterion to fix 
the minimum toi which ensures the validity of the Hol- 
stein model. Imposing the inequalities criterion AE^d > 
AE2D ^ AEiD we set the threshold value uji which turns 
out to be a function of the breathing mode energy wo and 
of the d— independent e-ph coupling go = g/Vd {g scales 
oc Vd): thus, at intermediate 5(0 (~ 1 — 1.5, in units of 
Wo) '^1 is — wo/2, while at strong go (> 2) cdi should be at 
least ~ 2a;o/3 in order to ensure the correct bandwidths 
trend. On the other hand, the intermolccular energies 
encounter the upper bound wi < luq given by the value 
of the coupling energy between the two atoms in the ba- 
sic unit of the molecular solid. Moreover, too large loi 
may invalidate strong coupling perturbative treatments 
of the Holstein model for three dimensional systems in 
fully adiabatic regimes [39]. 

With this caveat we study the polaron mass both in 
the Lang-Firsov and in the Modified Lang-Firsov scheme 
taking a lattice model in which first neighbors molecu- 
lar sites interact through a force constants pair poten- 
tial. Then, the d— dependent optical phonon spectrum 
is given by 

o ct ~\~ 'Y 1 / 

<^id{<i) = — ^ M 2ajcosq + 72 

^loici) = + ^x/a2+2a75(q)+7^(2 + Mq)) 

^loi^) = + ^V«^+2a7i(q)+7^(3 + Uq)) 

g{q) = cosqx + cosqy 
h{q) = 2cos{qx - qy) 
j(q) = cosqx + cosqy + cosqz 

l{q) = 2cos{qx - qy) + 2cos{qx - q^) + 2cos{qy - q^)) 

(0.14) 

where the intra-molccular force constant a and the 
inter-molecular first neighbors force constant 7 are re- 
lated to Wo and u)i by Wq = 2a/M and lo\ = "f/M re- 
spectively. M is the reduced molecular mass. In terms 
of Wo, the dimensionless parameter ztjiOQ defines the adi- 
abatic {zt/ujQ > 1) and the antiadiabatic {zt/ujQ < 1) 
regime. 

Second order perturbative theory introduces the po- 
laron mass m* dependence on the hopping integral t, 
hence on the adiabatic parameter, which would be absent 
in the first order Lang-Firsov theory. Generally, m* can 
vary with t/wo in two ways: m* becomes lighter either 
by increasing wo at fixed t or, by increasing t at fixed 
Wo. As the mass variation due to Wo is much stronger 
than that due to t, for a given adiabatic parameter, we 
may get difi'erent mass values according to the absolute 
values of wo and t. However, for sufficiently strong e-ph 
couplings which make the perturbative method applica- 
ble, the LF mass changes only slightly with t and sec- 
ond order corrections are small unless the intramolecular 
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phonon energies are low {llIq < ^OmeV) [39]. Hereafter 
we set ujQ = lOOmeV and select the adiabatic parameter 
by tuning t. 

In Figs.l, we plot the ratio of the one dimensional po- 
laron mass to the bare band mass against the e-ph cou- 
pling calculated both in the Lang-Firsov scheme and in 
the Modified Lang-Firsov expression. 

An intermediate regime 2t = loq is assumed in Fig. 1(a) 
while the intermolecular energy spans a range of weak to 
strong values. The striking different behavior between 
the LF and the MLF mass occurs for intermediate g 
while at very strong couplings the MLF plots converge, 
as expected, towards the LF predictions. The LF method 
overestimates the polaron mass ior g G [~ 1 — 2] (accord- 
ing to the value of wi) and mostly, it does not capture the 
rapid mass increase found instead in the MLF descrip- 
tion. Note that, around the crossover, the MLF polaron 
mass is of order ten times the bare band mass in the case 
oji = GOmeV. Large intermolecular energies enhance the 
phonon spectrum thus reducing the effective masses in 
both figures. In the MLF method, large uii tend also to 
smooth the mass behavior in the crossover region. 

Going to a fully adiabatic regime (see Fig. 1(b)) the 
discrepancies between LF and MLF plots are even more 
pronounced and the range of e-ph couplings in which 
the two methods converge shrinks considerably. There 
is scarce renormalization in the MLF curves up to the 
crossover which is clearly signalled by a sudden although 
continuous mass enhancement whose abruptness is signif- 
icantly smoothed for the largest values of intermolecular 
energies. In the antiadiabatic case shown in Fig. 1(c), 
the picture changes drastically and we recover a nearly 
coincident mass behavior in the LF and MLF methods 
throughout the whole range of couplings. The conver- 
gence is favoured at large lui . As above mentioned the LF 
plots show a strong resemblance in going from Fig. 1(a) to 
Fig. 1(c): infact, the LF method slightly depends on the 
hopping integral in ID systems with large intramolecular 
energy. The results we have displayed so far induce to 
reconsider the concept of self-trapping traditionally in- 
dicating an abrupt, but continuous, transition between 
an infinite size states at weak e-ph couplings and a finite 
(small) size polaron at strong e-ph couplings. According 
to the adiabatic polaron theory [40,11] there is no self- 
trapping event in one dimension as the polaron solution is 
always the ground state of the system. Instead, in higher 
dimensionality a minimum coupling strength is required 
to form finite size polarons, hence self-trapped polarons 
can exist at couplings larger than that minimum. As a 
shrinking of the polaron size yields a weight increase, the 
polaron mass behavior is accepted to be the most reli- 
able indicator of the self-trapping transition. The latter 
appears to us as a crossover essentially dependent on the 
degree of adiabaticity of the system and crucially shaped 
by the internal structure of the phonon cloud which we 
have modelled by tuning the intermolecular forces. We 
are then led to relocate the self trapping event in the pa- 
rameter space of ID systems admitting that also finite 
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size polarons can self-trap if a sudden change in their ef- 
fective mass occurs for some values of the e-ph couplings 
in some portions of the adiabatic regime. As fluctua- 
tions in the lattice distortions around the electron site 
are included in our variational wave function discontinu- 
ities in the polaron mass should not appear at the on- 
set of the transition [41]. Mathematically we select the 
crossover points through the simultaneous occurence of 
a maximum in the first logarithmic derivative and a zero 
in the second logarithmic derivative of the MLF polaron 
mass with respect to the coupling parameter: such in- 
flection points, corresponding to the points of most rapid 
increase for m* , are reported on in Figures 2, where the 
mass ratios are plotted for a wide choice of antiadiabatic 
to adiabatic regimes and a sizeable value of both in 
one, two and three dimensions. 

Some well known features of the antiadiabatic polaron 
landscape are confirmed by our analytical variational 
model in all dimensionalities: i) antiadiabatic polarons 
arc generally heavier than adiabatic ones although, at 
very strong couplings, the mass values converge at the 
Lang-Firsov results and ii) there is no self-trapping in 
the fully antiadiabatic regime as the electron and the 
dragged phonon cloud form a compact unit, a small po- 
laron, also at intermediate e-ph couplings. Then, the 
mass increase is smooth in the antiadiabatic regime. In- 
stead, in the more controversial [42] antiadiabatic to adia- 
batic transition region wc start to detect the signatures of 
the crossover which persist in the fully adiabatic regime 
and form a line of self-trapping events whose features 
however change considerably versus dimensionality. In 
ID (Fig. 2(a)), the crossover occurs for g values between 
^ 1.8 — 2.3 and the corresponding self trapped masses 
are of order ^5 — 50 times the bare band mass thus 
suggesting that relatively light small polarons can ex- 
ist in ID molecular solids with high phonon spectrum. 
The self-trapped mass values grow versus g by increasing 
the degree of adiabaticity and the incipience of the self- 
trapping line is set at the intermediate value 2t/ujQ = 1. 
We note that these findings are in good qualitative and 
quantitative agreement with refined variational results 
supporting the existence of self-trapped polarons also in 
ID. Although in the deep adiabatic regime we find a quasi 
step- like increase, the ID polaron mass is a continuous 
and derivable function of the e-ph coupling [23] . 

The two dimensional lattice introduces some signifi- 
cant novelties in the MLF mass behavior as shown in 
Fig. 2(b): i) at a given e-ph coupling and adiabaticity ra- 
tio, the 2D mass is lighter than the ID mass and the 
2D LF limit is attained at a value which is roughly one 
order of magnitude smaller than in ID; ii) the crossover 
region is shifted upwards along the g axis with the self 
trapping events taking place in the range, g ~ 2.2 — 2.6 
and the corresponding masses are of order ~ 5 — 10 times 
the bare band mass; iii) the curve connecting the self 
trapping points is parabolic with an extended descending 
branch starting at the intermediate value At/ijjQ= 1; iv) 
in the deep adiabatic regime, the lattice dimensionality 
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smoothens the mass increase versus g. The latter effect 
is even more evident in 3D, see Fig. 2(c), as there are no 
signs of abrupt mass increase even for the largest values 
of the adiabatic parameter. At the crossover, 3D masses 
are of order ~ 5 — 10 times the bare band mass with the 
self trapping points lying in the range, g ^ 2.5 — 2.9. At 
large couplings the effective mass over bare band mass 
ratio becomes independent of the t value and converges 
towards the LF vahie. In this region (and for the choice 
u)i = QOmeV) the 3D Lang-Firsov mass is one order 
of magnitude smaller than the 2D mass. As the coor- 
dination number grows versus dimensionality, large in- 
termolecular forces are more effective in hardening the 
3D phonon spectrum thus leading to lighter 3D polaron 
masses than 2D ones. 

In Figures 3 we plot the correlation functions Xo, Xi 
and X2 in ID (a) and 2D (b) respectively, as obtained 
by (13) for the adiabatic regime zt/uio = 2 with wo = 
lOOmeV. Two values, toi = AOmeV and wi = SOmeV, 
have been chosen to point out the role of the intermolec- 
ular forces in the transition between a large polaron at 
weak couplings and a small polaron at strong couplings. 
For sufficiently strong g values the LF limit is obtained, 
i.e. xa becomes 1 while xi and X2 become zero implying 
that the resulting polaron is an on site small polaron. 
The small to large polaron cross-over is manifested by a 
strong reduction of xo alongwith an enhancement in the 
values of xi and X2- By increasing wi, the crossover is 
slightly smoothed and shifted upwards along the g axis. 
Accordingly, xi and (to a lesser extent) X2 acquire some 
weight throughout a larger portion of e-ph coupling val- 
ues. As a main feature we note that the crossovers indi- 
cated by the correlation functions of the one dimensional 
system, for the two selected cases, occur at g/ujQ ^ 2 
and ~ 2.35 respectively. These values match the corre- 
sponding crossover points extracted by the polaron mass 
slopes. In two dimensions, the self-trapping transition 
takes place at larger (than in ID) g values and non local 
e-ph correlations persist in the adiabatic polaron up to 
g/uo ^ 3. The crossover is generally smooth and the 
softening effect of the intermolecular forces is more pro- 
nounced than in the one dimensional system. 

IV. CONCLUSIONS 

We have developed a variational analytical method to 

study the Holstcin polaron problem versus dimension- 
ality in the entire range of (anti)adiabatic parameters 
characterizing the molecular system. The essential role 
of the phonon dispersion in the Holstcin model has been 
accoimted for including the intermolecular interactions 
by means of a force constant approach. Unlike the tra- 
ditional Lang-Firsov scheme the Modified Lang-Firsov 
method permits to describe the fact that, in the inter- 
mediate and adiabatic regimes, the lattice deformation 
does not follow instantaneously the electron motion thus 
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leading to a spreading in the quasiparticle size. Under 
these circumstances we have examined the behavior of 
the polaron mass as a function of the strength of the e- 
ph coupHng and critically analysed the occurence of the 
self-trapping event signalling a shrinking of the polaron 
size in the real space. This crossover has been also mon- 
itored through the computation of the static e-ph cor- 
relation functions which provide a complementary tool 
corroborating our conclusions. Varying the adiabatic pa- 
rameter and selecting the points of most rapid increase 
for the effective mass we have found a set of self-trapping 
points originating, in ID and 2D, in the intermediate 
regime (zt/ujQ — 1) and continuing in the fully adia- 
batic regime. In 3D, the self-trapping events occur at 
zt/ too > 1. While, in one dimension, the curve connecting 
the inflection points in the adiabatic regime is a mono- 
tonic growing function of the e-ph coupling, in two and 
three dimensions we find distinctive parabola-like curves 
whose minima (of order ~ 5 times the bare band mass) 
arc located at larger g in higher d. Hence small polaron 
formation is favoured in low d whereas very large e-ph 
couplings are required to shrink the size of adiabatic po- 
larons in 3D. As intermolecular forces play a stronger 
role in more closely packed structures, lattice dimension- 
ality is expected to shape the polaron behavior. Infact, 
our results show that the crossover from large to small 
polarons is, in 2D and even more in 3D, smoother than 
in the case of the ID adiabatic polaron at a fixed value 
of intermolecular energy. Pointed out the quantitative 
differences in the polaron mass according to the dimen- 
sionality one should however notice a qualitative simi- 
larity in all dimensions regarding the occurence of the 
self-trapping event. Finally we observe that, although 
polaron masses become generally lighter in higher d, also 
in ID the effective mass over bare band mass ratio is 
^ 5 at the crossover when phonons and electrons com- 
pete on the energy scale. Small polarons having mobility 
properties may be therefore expected in low dimensional 
molecular systems with sufSciently strong intermolecular 
forces. 
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FIG. 1. Ratio of the one dimensional polaron mass to 
the bare band mass versus e-ph coupling according to the 
Lang-Firsov and the Modified Lang-Firsov methods. The 
adiabatic parameter is set at: (a) the intermediate value, 
2t/u)o = 1; (b) a fully adiabatic regime, 2t/Lj{) = 2; (c) an 
antiadiabatic regime, 2t/u}o = 0.25. u)o = lOOmeV and wi 
(in units meV) are the intramolecular and intermolecular en- 
ergies of the diatomic molecular chain respectively. 
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FIG. 2. Ratio of the Modified Lang-Firsov polaron mass to 
the bare band mass versus e-ph couphng in (a) ID, (b) 2D and 
(c) 3D. A set of twelve zt/uio values ranging from the antiadia- 
batic to the adiabatic regime is considered. From left to right: 
zt/ujo = 0.25, 0.5, 0.75, 1.0, 1.25, 1.5, 1.75, 2.0, 2.25, 2.5, 2.75, 3.0. 
ojo = lOOmeV^. The diamonds mark the occurence of the 
self-trapping event. 
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FIG. 3. (a) One dimensioned and (b) two dimensional static 
correlation functions versus e-ph coupling in the adiabatic 
regime, zt/ujo = 2. ujo = WOmeV. Two values of inter- 
molecular energies, uji = 40meV and uji = SOmeV, have 
been taken. 
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